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The hadron mass spectrum is calculated in lattice QCD using a novel fat-link clover fermion action in which only 
the irrelevant operators of the fermion action are constructed using smeared links. The simulations are performed 
on a 16^ X 32 lattice with a lattice spacing of a = 0.125 fm. We compare actions with n = 4 and 12 smearing 
sweeps with a smearing fraction of 0.7. The n = 4 Fat Link Irrelevant Clover (FLIC) action provides scaling 
which is superior to mean-field improvement, and offers advantages over nonperturbative 0{a) improvement. 



1. INTRODUCTION 

Understanding the generation of hadron mass 
from first principles has proved to be challeng- 
ing. The only method for deriving hadron masses 
directly from QCD is a numerical calculation on 
the lattice. The high computational cost required 
to perform accurate lattice calculations at small 
lattice spacings, however, has led to an increased 
interest in quark action improvement. To avoid 
the famous doubling problem, Wilson fermions [Q 
introduce additional terms which explicitly break 
chiral symmetry at 0{a). To extrapolate reliably 
to the continuum, simulations must be performed 
on fine lattices, which are very computationally 
expensive. The scaling properties of the Wilson 
action at finite a can be improved by introducing 
any number of irrelevant operators of increasing 
dimension which vanish in the continuum limit. 

The Sheikholeslami-Wohlert (clover) action 
introduces an additional irrelevant dimension-five 
operator to the standard Wilson [Q quark action, 

S'sw = Sw- ^"'^^^^ ■4>{x)af,„Ff,„ip{x) , (1) 

where S-w is the standard Wilson action and Csw 
is the clover coefficient which can be tuned to re- 
move 0{a) artifacts. Nonperturbative (NP) 0{a) 
improvement ||^ tunes Csw to all powers in 
and displays excellent scaling, as shown by Ed- 
wards et al. B, who studied the scaling prop- 



erties of the nucleon and vector meson masses 
for various lattice spacings (see also Section ^ be- 
low). In particular, the linear behavior of the NP- 
improved clover actions, when plotted against a^, 
demonstrates that 0{a) errors are removed. It 
was also found in Ref. |Q that a linear extrapo- 
lation of the mean- field improved data fails, indi- 
cating that 0{a) errors are still present. 

A drawback to the clover action, however, is the 
associated problem of exceptional configurations, 
where the quark propagator encounters singular 
behavior as the quark mass becomes small. In 
practice, this prevents the use of coarse lattices 
(/3 < 5.7 ~ a > 0.18 fm) [|j6|. Furthermore, the 
plaquette version of which is commonly used 
in Eq. (0), has large 0{a^) errors, which can lead 
to errors of the order of 10% in the topological 
charge even on very smooth configurations [[7| . 

The idea of using fat links in fermion actions 
was first explored by the MIT group ||^ and more 
recently has been studied by DeGrand et al. |^,|| , 
who showed that the exceptional configuration 
problem can be overcome by using a fat-link (FL) 
clover action. Moreover, the renormalization of 
the coefficients of action improvement terms are 
small. A drawback to conventional fat-link tech- 
niques, however, is that in smearing the links 
gluon interactions are removed at the scale of the 
cutoff. While this has some tremendous benefits, 
the short-distance quark interactions are lost. As 
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a result decay constants, which are sensitive to 
the wave function at the origin, are suppressed. 

A solution to these problems is to work with 
two sets of links in the fermion action. In the rel- 
evant dimension-four operators, one works with 
the untouched links generated via Monte Carlo 
methods, while the smeared fat links are intro- 
duced only in the higher dimension irrelevant op- 
erators. In this way the continuum limit of the 
theory is perfectly well defined. 

In this paper we present the first results of 
simulations of the spectrum of light mesons and 
baryons using this variation on the clover action. 
In particular, we will start with the standard 
clover action and replace the links in the irrel- 
evant operators with APE smeared |]l^, or fat 
links. We shall refer to this action as the Fat- 
Link Irrelevant Clover (FLIC) action. Although 
the idea of using fat links only in the irrelevant op- 
erators of the fermion action was developed here 
independently, suggestions have appeared previ- 
ously |l|]. To the best of our knowledge, this is 
the first report of lattice QCD calculations using 
this novel fermion action. 

2. GAUGE ACTION 

The simulations are performed using a tree- 
level C'(a^)-Symanzik-improved gauge action 
on a 16'^ X 32 lattice at /? = 4.60, providing a lat- 
tice spacing a = 0.125(2) fm determined from the 
string tension with i/a = 440 MeV. A total of 
50 configurations are used in this analysis, and 
the error analysis is performed by a third-order, 
single-elimination jackknife, with the per de- 
gree of freedom (A^df) obtained via covariance 
matrix fits. 

3. FAT-LINK IRRELEVANT FERMION 
ACTION 

Fat links HJ^] are created by averaging or 
smearing links on the lattice with their near- 
est neighbors in a gauge covariant manner (APE 
smearing). The smearing procedure [l^ replaces 
a link, C/^(x), with a sum of the link and a times 
its staples 
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+Ul{x — i'a)Ufj^{x — va)U,y{x — va + fj,a) 



followed by projection back to SU(3). We select 
the unitary matrix U^^^ which maximizes 

7^etr(C/;Lc/n) 

by iterating over the three diagonal SU(2) sub- 
groups of SU(3). We repeat this procedure of 
smearing followed immediately by projection n 
times. We create our fat links by setting a — 0.7 
and comparing n — A and 12 smearing sweeps. 
The mean-field improved FLIC action now be- 
comes 



FL 



Qt<Li rt 



z[Uq ) 



where F^^, is constructed using fat links, and 
where the mean-field improved Fat-Link Irrele- 
vant Wilson action is 
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with K = 1/(2™ + 8r). We take the standard 
value r = 1. The 7-matrices are hermitian and 

As reported in Table |l|, the mean-field improve- 
ment parameter for the fat links is very close to 
1. Hence, the mean-field improved coefficient for 
Csw is expected to be adequate^. In addition, 
actions with many irrelevant operators {e.g. the 
D234 action) can now be handled with confidence 

^Our experience with topological charge operators sug- 
gests that it is advantageous to include uq factors, even as 
they approach 1. 
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Table 1 

The value of the mean link for different numbers 
of smearing sweeps, n. 



n 


"0 







0.88894473 


0.62445197 


4 


0.99658530 


0.98641100 


12 


0.99927343 


0.99709689 



as tree-level knowledge of the improvement coeffi- 
cients should be sufficient. Another advantage is 
that one can now use highly improved definitions 
of F^jy (involving terms up to mJ^), which give 
impressive near-integer results for the topological 
charge 

In particular, we employ an ©(a**) improved 
definition of jl^ in which the standard clover- 
sum of four 1x1 Wilson loops lying in the fi, v 
plane is combined with 2x2 and 3x3 Wilson 
loop clovers. 

Work by DeForcrand et al. suggests that 
7 cooling sweeps are required to approach topo- 
logical charge within 1% of integer value. This 
is approximately 16 APE smearing sweeps at 
a = 0.7 ||l^. However, achieving integer topo- 
logical charge is not necessary for the purposes 
of studying hadron masses, as has been well es- 
tablished. To reach integer topological charge, 
even with improved definitions of the topologi- 
cal charge operator, requires significant smooth- 
ing and associated loss of short-distance informa- 
tion. Instead, we regard this as an upper limit on 
the number of smearing sweeps. 

Using unimproved gauge fields and an unim- 
proved topological charge operator. Bonnet et 
al. found that the topological charge settles 
down after about 10 sweeps of APE smearing at 
a = 0.7. Consequently, we create fat links with 
APE smearing parameters n = 12 and a = 0.7. 
This corresponds to ^ 2.5 times the smearing 
used in Refs. |^,^. Further investigation reveals 
that improved gauge fields with a small lattice 
spacing (a — 0.125 fm) are smooth after only 
4 sweeps. Hence, we perform calculations with 



4 sweeps of smearing at a — 0.7 and consider 
n = 12 as a second reference. Table |l| lists the 
values of Uq^ for n = 0, 4 and 12 smearing sweeps. 

A fixed boundary condition is used for the 
fermions by setting 

Ut{x, nt) = and Uj^{x, nt) = M x (5) 

in the hopping terms of the fermion action. The 
fermion source is centered at the space-time lo- 
cation {x,y,z,t) = (1,1,1,3), which allows for 
two steps backward in time without loss of sig- 
nal. Gauge-invariant gaussian smearing p^ ] in 
the spatial dimensions is applied at the source to 
increase the overlap of the interpolating operators 
with the ground states. 

4. RESULTS 

Hadron masses are extracted from the Eu- 
clidean time dependence of the calculated two- 
point correlation functions. The effective masses 
are given by 

M{t + 1/2) = log[G(i)] - \og[G{t + 1)] . (6) 

The critical value of k, Kc, is determined by lin- 
early extrapolating as a function of to 
zero. We used five values of quark mass and the 
strange quark mass was taken to be the second 
heaviest quark mass. 

Effective masses (^ are calculated as a func- 
tion of time and various time-fitting intervals 
are tested with a covariance matrix to obtain 
X^/A^DF- Good values of x^/A^df are obtained 
for many different time-fitting intervals as long 
as one fits after time slice 8. All fits for this ac- 
tion ("FLIC4") are therefore performed on time 
slices 9 through 14. For the Wilson action and the 
FLIC action with n = 12 ("FLIC12") the fitting 
regimes used are 9-13 and 9-14, respectively. 

The behavior of the p, nucleon and A masses as 
a function of squared pion mass is shown in Fig. || 
for the various actions. The first feature to note is 
the excellent agreement between the FLIC4 and 
FLIC12 actions. On the other hand, the Wilson 
action appears to lie somewhat low in compar- 
ison. It is also reassuring that all actions give 
the correct mass ordering in the spectrum. The 
value of the squared pion mass at m-,^lmp ~ 0.7 
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Figure 1. Masses of the nucleon, A and p meson 
versus for the FLIC4, FLIC12 and Wilson 
actions. 
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is plotted on the abscissa for the three actions as 
a reference point. This point is chosen in order to 
allow comparison of different results by interpo- 
lating them to a common value of m^/mp = 0.7, 
rather than extrapolating them to smaller quark 
masses, which is subject to larger systematic and 
statistical uncertainties. 

The scaling behavior of the different actions 
is illustrated in Fig. ^. The present results for 
the Wilson action agree with those of Ref. 
The first feature to observe in Fig. ^ is that ac- 
tions with fat-link irrelevant operators perform 
extremely well. For both the vector meson and 
the nucleon, the FLIC4 action performs system- 
atically better than the FLIC12. This suggests 
that 12 smearing sweeps removes too much short- 
distance information from the gauge-field configu- 
rations. On the other hand, 4 sweeps of smearing 
combined with our O(a^) improved F^i, provides 
excellent results, without the fine tuning of Csw 
in the NP improvement program. Notice that for 
the p meson, a linear extrapolation of the mean- 
field improved clover points in Fig. ^ indicates 
that there is better improvement when using fat 
links in the irrelevant operators. While there are 
no NP-improved clover plus improved glue simu- 
lation results at a^a ^ 0.08, the simulation results 



Figure 2. Nucleon and vector meson masses for 
the Wilson, NP-improved and FLIC actions ob- 
tained by interpolating our results of Fig. to 
m-^jmp = 0.7. Results from the present simula- 
tions are indicated by the solid points. The fat 
links arc constructed with n = 4 (solid squares) 
and n = 12 (stars) smearing sweeps at a = 0.7. 



that are available indicate that the fat-link results 
also compete well with those obtained with a NP- 
improved clover fermion action. 

Finally, we compare the convergence rates of 
the different actions by comparing the number of 
Stabihsed Biconjugate Gradient (BiCGStab) Q 
iterations required to invert the fermion matrix. 
The FLIC 12 action converges at a slightly faster 
rate to the Wilson fermion action and the FLIC4 
action is the clear winner due to the small num- 
ber of BiCGStab iterations required to invert the 
fermion matrix. This provides great promise for 
performing simulations at quark masses closer to 
the physical values. 
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5. CONCLUSIONS 

We have examined the hadron mass spec- 
trum using a novel Fat Link Irrelevant Clover 
(FLIC) fermion action, in which only the irrele- 
vant, higher-dimension operators involve smeared 
links. One of the main conclusions of this work is 
that the use of fat links in the irrelevant opera- 
tors provides excellent results. Fat links promise 
improved scaling behavior over mean-field im- 
provement. This technique also solves a sig- 
nificant problem with 0{a) nonperturbative im- 
provement on mean field-improved gluon config- 
urations. Simulations are possible and the results 
are competitive with nonperturbative-improved 
clover results on plaquette-action gluon configu- 
rations. We have found that minimal smearing 
holds the promise of better scaling behavior. Our 
results suggest that too much smearing removes 
relevant information from the gauge fields, lead- 
ing to poorer performance. Fermion matrix inver- 
sion for FLIC4 is more efficient and results show 
no sign of exceptional configuration problems. 

This work paves the way for promising future 
studies. It will be of great interest to consider dif- 
ferent lattice spacings to further test the scaling of 
the fat-link actions. Furthermore, the exceptional 
configuration issue can be explored by pushing 
the quark mass down to lower values. The n = 4 
FLIC action holds great promise for circumvent- 
ing this issue as evidenced by the relative ease 
with which one can invert the fermion matrix. A 
study of the spectrum of excited hadrons using 
the fat-link clover actions is currently in progress 
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